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Radio-frequency electric-dipole transitions between nearly degenerate, opposite 
parity levels of atomic dysprosium (Dy) were monitored over an eight-month 
period to search for a variation in the fine-structure constant, a. The data 
provide a rate of fractional temporal variation of a of (— 2.4it2.3) X 10~^® yr~^ 
or a value of (—7.8 it 5.9) X 10~^ for ka, the variation coefficient for o in a 
changing gravitational potential. All results indicate the absence of significant 
variation at the present level of sensitivity. We also present initial results on 
laser cooling of an atomic beam of dysprosium. 
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1. Introduction 

A component of Einstein's equivalence principle (EEP) is local position in- 
variance (LPI), which states that the laws of physics, including the values 
of fundamental constants, should be independent of space and time. Mod- 
ern theories attempting to unify gravitation with the other fundamental 
interactions allow, or even predict, violations of EEP,^ which has sparked 
searches for violation of LPI, and hence EEP, through searches for temporal 
and spatial variation of fundamental constants. 

Various studies have reported results of searches for a temporal vari- 
ation of the fine-structure constant (a = e^/hc) over cosmological time 
scales of 10^° years, geological time scales of 10^ years, and present 
day laboratory searches over the course of years. ^"'^^ In contrast to stud- 
ies involving analyses of processes that have occurred billions of years ago, 
the results from laboratory searches are easier to interpret since the ex- 
periments are repeatable, and systematic uncertaintities can be studied by 
changing experimental conditions. The best limit on a temporal variation 
that is independent of assumptions regarding other constants (published 
after our main result in Ref. 14) was obtained by monitoring the ratio of 
A1+ and Hg+ optical transition frequencies.^^ 

Another type of search for an LPI violation is a "null" gravitational 
red-shift experiment where two clocks with different composition are com- 
pared side by side in a changing gravitational potential. Laboratory clock 
comparisons can be used for this type of test due to the eccentricity of 
Earth's orbit around the Sun, which leads to a small oscillatory component 
of the gravitational potential with a period of a year. In earlier work,^^'^''"^^ 
clocks of different types were compared, and the ratios of the clock rates 
were analyzed for a possible correlation with the gravitational potential, 
which led to bounds on parameters that characterize structure-dependent 
modifications to the clock frequencies. At the time of publication, we were 
able to present the best limit on fc^, the linear- variation coefficient for a in a 
changing gravitational potential. ^° A better result has since been published 
in Ref. 21. 

Advances in direct comparison of single-ion optical or neutral optical 

lattice clocks through frequency-comb metrology promise significant im- 
provements in sensitivity by up to three orders of magnitude, as well as 
simplification in interpretation of the results since such comparisons di- 
rectly probe a-variation, independent of other fundamental constants (see, 
for example Ref. 15). In this article, we present results obtained with an 
alternative method of competitive sensitivity, which is also independent of 
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Fig. 1. Relevant levels and transitions in atomic dysprosium. Population and detection 
scheme: Level B is populated in a three-step process. The atoms are excited to level / in 
the first two steps using 833- and 699- nm laser light (solid arrows). The third step is a 
spontaneous decay (labeled by a short-dashed arrow) from level / to level B with ~ 30 % 
branching ratio. The population is transferred to level A by the rf electric field (wavy 
line). Atoms in level A decay to level c and then to the ground state. The fluorescence 
from the second decay (564 nm) is used for detection. Inset on the right: Hyperfine 
structure of levels A and B for ^"^^Dy. Zero energy is arbitrarily chosen to coincide with 
the lowest hyperfine component. 

other fundamental eonstants. This method utihzes rf electric-dipole transi- 
tions between nearly degenerate electronic levels in atomic Dy. 

2. Effect of a-variation in dysprosium 

Tests for variation of a in atomic systems rely upon the fact that the sensi- 
tivity of energy levels to a is different for different energy levels. The total 
energy of a level can be written as 

E^hi^ = Eo + q{a^/al-l), (1) 

where Eq is the present-day energy, ao is the present-day value of the 
fine-structure constant, and g is a coefficient which determines the sen- 
sitivity to variation of a. A calculation,^^ utilizing relativistic Hartree- 
Fock and configuration interaction methods, found values of q for two 
nearly degenerate opposite-parity levels in atomic dysprosium (Fig. 1, lev- 
els A and B) that are large and of opposite sign. For the even-parity level 
(level A), QA/hc ~ 8 x 10'^ cm^"'^, while for the odd-parity level (level B), 
qs/hc 2± -25 X fO^ cm"^ 
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Fig. 2. Experimental setup (not shown to scale): a) Atomic beam produced by eflfusive 
oven source at ~ 1500 K; b) atomic-beam collimator; c) oven light collimator; d) cylin- 
drical lenses to diverge laser beams; c) interaction region of atoms with the electric field 

enclosed in a magnetic shield (not shown); f) cylindrical mirror to collect fluorescent 
light; g) Incite light pipe; h) interference filter; and i) short-pass filter. 



3. Experimental Technique 

A unique aspect of the Dy system is that the directly observable quantity 
is the energy difference of the levels, vb—i^a, due to the fact that an electric- 
dipole transition can be induced between them. The time variation of the 
transition frequency between levels A and B is related to temporal variation 
of a by 



Au = ub - = 2 , (2) 

ha 

where \2{qB - qA)/h\ ~ 2.0 x 10^^ Hz.^^ For instance, \a/a\ = IQ-'^^/yr 
implies |Az/| ~ 2 Hz/yr. 

A preliminary analysis of statistical and systematic uncertainties shows 
that the measurement of the transition frequency and the control of the 

systematics is feasible at a mHz level, which corresponds to an ultimate 
sensitivity of \a/a\ ~ 10~^*/yr for two measurements separated by a year's 
time.^^ A mHz resolution on a transition frequency of 1 GHz requires a 
relatively modest fractional stability of 10~^^ for the reference frequency 
standard. This also means that the results are insensitive to variation of 
the Cs reference frequency due to changes in the values of fundamental 
constants, which have experimental upper limits~ 10~^^ yr~^ (see, for ex- 
ample, Ref. 10). 

The transitions utilized in the experiment are the 235-MHz (J = 10 
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Fig. 3. Measured transition frequencies for (a) the 235-MHz and (b) 3.1-MHz transitions 
over an eight-month period. The data have been corrected for coUisional shifts. The sohd 
lines are the least-squares linear fits to the data. The dashed lines are the least-squares 
fits to a constant function. The apparent correlation between the transition frequencies 
for the two isotopes are due to the stray magnetic field effect. 

^ J = 10) transition of the isotope ^^^Dy and the 3.1-MHz (F = 10.5 
^ F = 10.5) transition of ^'^^Dy. The long-hved level B (hfetime^'* r > 
200 /is) is populated via three transitions. Lasers at 833 nm and 669 nm 
are used for the first two transitions to a high-lying upper state, which then 
spontaneously decays to level B as shown in Fig. 1. The transition between 
levels A and B is induced by an applied rf electric field. The atoms decay 
from level A in two steps, and the fluorescence at 564 nm is monitored to 
detect the rf transition. 

The rf-generation and detection system is discussed in Ref. 25. The rf 
field is frequency modulated by the 10 kHz reference output of a lock-in 
amplifier with a modulation index of 1. In order to reduce lineshape asym- 
metries due to drifts, the ratio of the first- and second-harmonic outputs of 
the lock-in amplifier is used to measure the transition frequency, which is 
extracted from the ratio by a two-step process described in detail in Ref. 25. 

4. Results and Analysis 
4.1. Temporal Variation 

We have measured the two rf transition frequencies in the course of seven 
runs over eight months. Figure 3 shows the results of these measurements 
corrected for coUisional shifts. The present uncertainty is dominated by 
systematic effects primarily due to polarization imperfections coupled to 
the residual magnetic field, coUisional shifts, and rf electric-field inhomo- 
geneities. A detailed analysis of these systematic effects and a discussion of 
feasible levels of stability are given in Ref. 23. Error bars are due both to sta- 
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Fig. 4. a) The true anomaly, 0, is the angle subtended from perihelion, the point of 
closest approach. The eccentric anomaly, e, is the angle between perihelion and the 
position of the Earth in its orbit projected onto the auxiliary circle of the ellipse (the 
eccentricity of the ellipse is exaggerated for clarity), b) The change in gravitational 
potential of the Sun at the Earth due to the ellipticity of the orbit. The points indicate 
the dates on which data were taken. 



tistical and systematic uncertainties which have been added in quadrature. 
A least-squares hnear fit to the data points gives slopes of —0.6 ± 6.5 Hz/yr 
and 9.0 ±6.7 Hz/yr for the 235-MHz and 3.1-MHz transitions, respectively. 
Using an analysis presented in Ref. 14, we arrive at the final result of 
aja — (— 2.4± 2.3) x 10^^^ yr~^, consistent (Icr) with no variation at the 
present level of sensitivity. This result is slightly improved over what is pre- 
sented in Ref. 14 due to recently improved calculations of the sensitivity 
coefficients, qA and gs.^^ 



4.2. Gravitational-Potential Dependence 

The data used in the previous analysis can also be analyzed for the depen- 
dence of a on variations of a gravitational potential. The data were taken 
over the course of several months, and because the Earth's orbit has a 
slight eccentricity, each measurement took place at a different gravitational 
potential. 

The fractional change in a as a function of gravitational potential, U{t), 
can be parameterized as^^ 

^ = t.^. (3) 

where ka is the variation coefficient for a in a changing gravitational poten- 
tial, c is the speed of light, and AU{t) is the oscillatory component of the 
sun's gravitational potential at the Earth (see Fig. 4 and Ref. 20). The mea- 
sured frequencies for each isotope are fitted to the gravitational potential 
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Fig. 5. Same data as in Fig. 3. Solid lines are least-squares fit to a function given by 
Eq.fefeq:fit where AU is the oscillatory component of the gravitational potential of the 
Sun at Earth with a period of a year. 

using a two-parameter least-squares fit given by 

6{Aiy) - i'^ = xo 5^+a;i, (4 

where v* is an arbitrary reference frequency, and xq and xi are the fit 
parameters. The parameter xi accounts for the offset due to the reference 
frequency while the parameter xq determines the correlation between the 
change in gravitational potential and the transition frequency. 

The parameters xq and xi were obtained by a least-squares fit to the 
data shown in Fig. 5. These values can be used to calculate the constraint 
on the parameter ka from Eq. (3). Substituting Eq. (2) into Eq. (3), we get 
a relation similar to Eq. (4), 

.,A.).(2i£^)..^, (5, 

The fit parameter yields ka = (-1.0 ± 8.2) x 10"^ for the 235-MHz 
transition, k^ = (-14.6 ±8.4) x 10"'^ for the 3.1-MHz transition. The final 
result of ka — (—7.8 ± 5.9) x 10^® is obtained from the analysis presented 
in Ref. 20. In addition, our value of ka, combined with other experimental 
constraints, is used to extract the first limits on kg and kq. These coefficients 
characterize the variation of nie/mp and mq/mp in a changing gravitational 
potential, where mg, rup, and rriq are electron, proton, and quark masses. 
The results are ke = (4.3 ±3.5) x lO'^ and fc, = (5.9 ±4.7) x 10"^ As was 
noted in Sec. 4.1, these results are improved over those presented in Ref. 20 
due to a more accurate calculation of sensitivity coefficients.^^ 



5. Laser Cooling of Dysprosium 

Recently, we have started to explore the possibility of laser cooling dyspro- 
sium, with the goal of generating a well coUimated beam for our experiment. 
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Fig. 6. Scans of the ^^^Dy peak in the 658-nin transition with (hollow points) and with- 
out (solid points) the 421-nrn pushing light incident on the atomic beam. In the scan on 
the left the 421-nm light was tuned to the center of the ^^'*Dy resonance, and in the scan 
on the right the 421-nm light was detuned from resonance by half the transverse doppler 
width. In both scans an alteration of the lineshape, and hence velocity distribution, is 
clearly apparent. 

This is expected to increase beam brightness, reduce sensitivity to spatial 
magnetic- and electric-field inhomogeneities, and allow us to characterize 
shifts associated with Dy-Dy collisions. 

While dysprosium has been magnetically trapped before, there is no 
previously known effort to perform laser cooling. After a review of literature 
concerning excited state lifetimes, -^^'^^ it was decided that the best option 
to attempt laser cooling was to use the 421.291-nm, 4/^''6s^ (J = 8) — > 
4/^°6s6p (J = 9) transition. With an upper state lifetime of 4.8 ns, this 
transition is the strongest recorded cycling transition in the visible spectrum 
of Dy 1. 

To generate the 421-nm light, 300 mW of 842-nm light from a 
Ti:Sapphire ring laser (Coherent 899) is frequency-doubled by a 1x2x10 
mm^ periodically poled potassium titanyl phosphate (PPKTP) crystal 
(Raicol Crystals Ltd.) located within a resonant bow-tie cavity. The max- 
imum output with this resonator-enhanced doubling scheme is approxi- 
mately 90 mW. 

In our initial attempt we choose to "push" the atomic beam, as cool- 
ing required more changes to the experimental set-up. To detect changes 

in the atomic beam, we use a parallel probe beam 1-cm downstream from 
the pushing beam. Both the pushing and probing laser beams intersect 
the atomic beam at ninety degrees. The 658.12-nm light for the probe 
laser is generated by a ring dye laser with DCM dye and resonant with 
the 4/i°6s2(J = 8) ^ 4:f5d6s'^{J = 7) transition. The linewidth of this 
transition^^ (~100 kHz) is much less than the Doppler width (40 MHz), 
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therefore the lincshapc of the probe fluorescence is dominated by the trans- 
verse velocity distribution of the atomic beam. A more detailed description 
of the atomic beam source is found in Ref. 30. 

The results of the experiment are shown in Fig. 6. For the first scan 
the 421-nm light is resonant with the center of the transverse velocity dis- 
tribution, resulting in a shift and broadening of the probe lineshape. In 
the second scan the 421-nm is red-detuned from the center of the veloc- 
ity distribution, resulting in an increase in amplitude and narrowing of the 
lineshape. In both scans the power in the pushing beam is approximately 60 
mW in a 1.1-cm diameter, coUimated beam. Work is currently underway to 
modify the experiment for the purposes of more efficient transverse cooling. 

Laser cooling of dysprosium is not trivial because of many available 
decay channels of the upper state of the cooling transition. However, our 
estimates (confirmed by more detailed atomic calculations^^) show that 
the overall branc;liing ratio into all excited states is < 10^"'. In future work, 
we will determine the branching ratios and devise repumping schemes as 
needed. In principle magneto-optical trapping, or even Bose-Einstein con- 
densation may be possible. The latter would be of interest to the degenerate 
atomic gas commrmity due to the large ground-state angular momentum 
and magnetic moment of Dy. 



6. Conclusion 

We have presented improved results of a measurement of time variation 
ofa and a correlation with a changing gravitational potential using atomic 
dysprosium. In our case the interpretation does not require comparisons 
with different measurements to eliminate dependence on other constants, 
and provides an alternative to measurements that utilize state-of-the-art 
atomic optical frequency clocks. In addition, the data have been used to 
extract a limit on the gravitational variation of a. We have also presented 
results on the first effort to laser cool dysprosium. 

A significant improvement in our results is expected from a newly op- 
eration apparatus, which provides better control over polarization imper- 
fections coupled to the residual magnetic field, collisional shifts, and rf 
electric-field inhomogeneities, as well as other systematic effects expected 
to be important to achieve better than 1-Hz sensitivity. Ultimately, mHz- 
level sensitivity may be achievable with this method. 
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